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ABSTRACT 

A carbon  nanotube  pulling  through  a polyethylene  matrix  was  simulated  using  molecular 
dynamics.  The  interfacial  sliding  was  characterized  in  terms  of  a nanoscale  friction  model,  which 
is  parametrized  from  the  molecular  dynamics  simulation,  and  involves  determining  the  critical 
pull-out  force  on  the  nanotube  and  the  effective  viscosity  at  the  nanotube/polymer  interface. 
Comparison  was  made  of  the  pull-out  behavior  of  non-bonded  and  functionalized  nanotube 
composites.  Chemical  bonds  between  the  polymer  and  the  nanotube  increased  the  critical  pull- 
out force,  the  resistance  to  interfacial  sliding,  and  the  interfacial  viscosity. 


INTRODUCTION 

In  polymer-nanotube  composites,  the  strength  of  the  interface  between  the  polymer  and  the 
nanotube  (NT)  may  be  enhanced  by  chemically  bonding  the  NTs  to  the  polymer.  Chemical 
functionalization  of  NTs  remains  a challenge  [1,  2],  but  in  some  polymer-NT  composites,  where 
significant  load  transfer  between  the  polymer  and  the  nanotube  has  been  observed,  chemical 
bonding  with  the  NT  may  occur  [3], 

In  traditional  fiber-reinforced  composites,  the  fiber-matrix  interfacial  shear  strength  is  often 
characterized  by  fiber  pull-out  experiments  [4],  Analogous  nanoscale  pull-out  experiments  have 
been  carried  out  with  atomic  force  microscopy  for  the  sliding  of  NTs  in  multi-walled  NTs  [5] 
and  in  bundles  of  single-walled  NTs  [6],  A simulated  pull-through  test  has  also  been  devised  for 
polymer-NT  composites  via  molecular  dynamics  (MD)  to  determine  the  interfacial  shear  strength 
while  a force  is  applied  to  the  NT.  These  simulations  have  indicated  an  increase  in  the  interfacial 
shear  strength  with  the  addition  of  chemical  bonds  between  the  polymer  and  a NT  [7], 
Furthermore,  the  simulations  have  been  extended  to  study  the  interfacial  sliding  of  the  NT  during 
pull-out  in  non-bonded  composites,  and  to  characterize  the  effective  interfacial  viscosity  due  to 
the  van  der  Waals  interactions  of  the  polymer  and  the  NT  [8,  9], 

In  this  work,  the  interfacial  sliding  during  pull-out  of  a NT  chemically  bonded  to  a crystalline 
polymer  matrix  was  modeled  with  MD  simulations.  The  objective  was  to  show  some  of  the 
mechanisms  which  may  serve  to  enhance  its  ability  to  transfer  mechanical  load  within  the 
composite.  First,  a description  of  the  composite  structure  will  be  given  with  details  of  the  MD 
simulations.  Then  stages  encountered  by  the  NT  during  interfacial  sliding  will  be  analyzed. 
Finally,  the  sliding  velocity  will  be  interpreted  by  using  the  interfacial  viscosity  model  [8, 9], 

MOLECULAR  DYNAMICS  SIMULATION 

The  simulated  polymer-NT  composite  consisted  of  a (10,10)  carbon  nanotube  embedded  in  a 
crystalline  polyethylene  (PE)  matrix  (Figure  la).  The  NT  was  chemically  bonded  at  three 
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locations  to  one  of  the  nearby  polymer  chains  (Figure  lb).  Each  attachment  included  a direct 
bond  to  the  NT,  two  CH2  groups  and  a direct  bond  to  the  chain  (Figure  Id).  In  the  MD 
simulations,  all  the  chemical  bonding,  including  the  chemical  bonds  between  the  NT  and  the  PE, 
was  represented  with  the  same  hydrocarbon  potential  [10].  The  non-bonded  interactions  were 
represented  with  the  Lennard-Jones  potential  [11],  The  system  was  replicated  across  periodic 
boundaries  in  each  of  the  3-dimensions,  making  the  NT  and  the  PE  chains  infinitely  long.  This 
composite  was  similar  to  the  crystalline  polyethylene  used  in  previous  work  [7],  and  the  detailed 
chemical  description  is  kept  for  consistency. 

To  simulate  the  interfacial  sliding  of  the  NT,  force  was  added  to  each  atom  of  the  NT  in  the 
axial  direction.  As  the  system  was  periodic,  the  NT  moved  through  the  polymer  into  the  next 
periodic  image  when  it  began  to  slide.  The  applied  force  was  incremented  over  time.  The 
simulation  is  run  for  about  210,000  MD  steps  of  0.05  fs  each.  The  size  of  the  force  increment 
was  increased  by  a factor  of  2.5  at  105,000  steps.  The  displacements  and  velocities  of  the  NT, 
the  attached  chains,  and  the  nearby  chains  were  monitored  under  increasing  force.  The  attached 
chain  exhibited  a breakage  before  the  initial  force  loading  due  to  reactivity  included  in  the 
hydrocarbon  potential,  and  retained  that  configuration  throughout  the  rest  of  the  simulation 
(Figure  lc).  The  initial  configuration  also  included  attachments  to  a second  chain  on  the  opposite 
side  of  the  NT,  but  this  attached  chain  did  not  survive  past  40,000  steps.  Only  a third  of  this 
chain  survived  the  initial  NT  pull-out  force  intact.  This  incident  had  marginal  effect  of  the 
investigated  interfacial  sliding  after  the  onset  of  pull-out;  that  sliding  was  dominated  by  the 
friction  along  the  entire  surface  area  and  the  chemical  bonding  of  the  first  chain. 


Figure  1 . The  structure  of  the  functionalized  polymcr-nanotubc  composite  (a)  cross-sectional 
view,  (b)  side  view  of  the  starting  structure  with  links  1 and  2 and  (c)  the  actual  structure  of 
NT  and  the  attached  chain  with  one  breakage  (d)  close-up  of  an  attachment. 


RESULTS  AND  DISCUSSION 
Interfacial  Sliding 

The  velocity  profile  of  the  NT  under  load  is  plotted  in  Figure  2.  Several  stages  were 
observed  in  the  interfacial  sliding  of  the  NT  through  this  composite.  However,  the  intcrfacial 
sliding  of  the  NT  in  this  composite  is  significantly  different  from  that  observed  previously  for  the 
non-bonded  PE/NT  composite  [8,  9],  Initially,  the  NT  experienced  thermal  vibrations  and 
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interacts  with  the  adjacent  polymer  chains  in  equilibrium  at  300  K,  At  0,17  nN  enough  force 
was  exerted  on  the  NT  for  interfacial  sliding  to  begin.  This  initial  pull-out  force  was  higher  than 
that  observed  for  the  onset  of  NT  sliding  in  a non-bonded  composite.  The  next  difference 
between  the  composites  occurred  when  the  NT  velocity  peaked  just  past  20,000  steps  and  then 
decreased  again  almost  to  zero.  At  this  point,  the  NT  was  under  the  constraint  of  the  chemical 
bonds  attaching  it  to  the  chain.  Increasing  the  force  on  the  NT  beyond  this  point  caused  it  to  drift 
about  0.5  A as  shown  by  the  NT  displacement  in  Figure  3.  At  about  57,000  steps,  the  NT  began 
to  accelerate  on  average  until  about  141,000  steps.  Although  the  size  of  the  force  increment  was 
increased  in  this  range,  the  relation  of  force  to  velocity  maintained  the  same  slope  (Figure  4). 
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GURE  2,  Velocity  of  the  center  of  mass  of  a carbon  NT  during  the  simulated  pull-, 
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Figure  3.  Displacement  of  individual  carbon  atoms  of  the  NT  during  the  simulated  pull-out. 


FIGURE  4.  The  force-velocity  relation  for  the  simulated  NT  pull-out. 


388 


In  this  range,  one  of  the  links  (Link  1 in  Figure  (lb))  of  the  chain  engaged  part  of  the  chain 
nearby  to  move  with  the  NT.  The  chain  was  loose  enough  that  the  part  near  Link  1 moved  with 
the  NT  while  the  rest  of  chain  progressively  slowed  down  in  the  vicinity  of  Link  2,  Link  2 was 
moving  with  the  NT,  but  it  was  loose  enough  that  it  did  not  engage  its  part  of  the  chain  into 
motion.  The  linkages  and  the  chains  had  significantly  more  degrees  of  freedom  and  could  move 
differently  relative  to  one  another  and  the  NT.  A transition  took  place  at  about  141,000  steps, 
where  Link  2 broke  releasing  its  part  of  the  chain.  This  transition  was  followed  by  two  peaks  in 
NT  velocity  (Figure  2)  indicating  continued  acceleration  of  the  NT  and  the  attached  chain. 

In  Figure  3,  the  displacement  of  individual  carbon  atoms  associated  with  the  sliding  NT,  the 
bonded  chain  and  a non-bonded  chain  are  shown  as  the  applied  force  increased  up  to  3 nN,  The 
bonded  polymer  chain  moved  along  with  the  NT  pulled  by  force,  while  a non-bonded  chain 
vibrated  around  its  initial  position  without  significant  motion  in  the  direction  of  NT  axis. 

Interfacial  viscosity 

The  NT  pull-out  process,  in  general,  and  interfacial  sliding,  in  particular,  can  be  described  by 
an  interfacial  friction  model  [9],  which  is  analogous  to  a Bingham-type  model  for  viscoplastic 
failure  of  the  interface  including  the  concept  of  interfacial  viscosity  [12,13],  The  model 
describes  the  dependence  of  the  average  applied  force,  on  the  average  NT  velocity,  <V:>, 

in  the  z-direction  of  the  NT  axis: 

a)  <f>P„n=fo+Zrff  <y:  >. 

where  fa  is  the  critical  pull-out  force  and  Zjj  is  the  effective  viscosity  coefficient,  which  can  be 
obtained  from  the  MD  simulation.  The  effective  viscosity,  of  the  NT/polymer  interface  is 
defined  as  a proportionality  constant  in  the  force-velocity  relation  that  characterizes  the 
interfacial  momentum  transfer.  NT  velocity  variations  shown  in  Figure  2 are  plotted  in  Figure  4 
for  every  incremental  force  loading.  Once  the  critical  pull-out  force, /o,  at  0.17  nN  was  reached, 
the  force-velocity  dependence  could  be  approximated  by  a linear  relation  (the  solid  line).  The 
inverse  of  the  slope  of  the  resulting  line  is  the  coefficient  %eg,  which  is  2.3+0.2  (nN  ps)/A,  This 
slope  can  be  estimated  from  dashed  lines  like  those  shown  in  Figure  4,  which  bound  the  NT 
velocity  variations  during  NT  sliding. 

The  average  applied  force,  in  equation  (1)  is  related  to  the  average  shear  stress,  <t>, 

for  sliding  interactions  by  the  force  balance:  <r>  » <f>n„n/A,,  where  A,  is  the  interfacial  area 
between  the  NT  and  polymer.  During  the  NT  pull-out,  the  applied  force  is  balanced  by  the 
interfacial  friction  characterized  by  the  effective  viscosity,  and  the  induced  shear  stress,  <t>. 
The  effective  viscosity  coefficient,  Xefft  is  related  to  the  effective  viscosity  via 

(2) 

A,. 

where  h.llw  is  the  van  der  Waals  separation  between  the  stationary  polymer  and  the  sliding  NT. 
Therefore,  the  effective  viscosity,  is  0.27  cP,  which  is  50%  higher  than  that  of  the  non- 
bonded  NT  composite,  0.18  cP  [8, 9], 
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CONCLUDING  REMARKS 


The  intcrfacial  sliding  of  a carbon  nanotubc  (NT)  chemically  bonded  to  a local  chain  in  a 
polymer  matrix  is  simulated  with  molecular  dynamics  (MD)  and  the  results  fit  an  interfacial 
friction  model  for  NT  pull-out.  The  interfacial  sliding  in  this  composite  differs  from  a non- 
bonded  composite  in  at  least  two  major  ways.  First,  the  NT  is  observed  to  be  constrained  by  the 
chemical  linkages.  Second,  as  the  linkages  are  flexible,  they  engage  separately  in  motion  with 
the  NT  at  different  times  and  force  loadings.  These  differences  result  in  higher  interfacial  shear 
strengths  and  higher  interfacial  viscosities  than  in  the  non-bonded  composite,  yielding  qualitative 
information  for  the  design  of  NT-reinforced  polymer  composites. 
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